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A kinetic theory for radiation interacting with sound waves in an ultrarelativistic 
electron-positron plasma is developed. It is shown that the effect of a spatial spectral 
broadening of the electromagnetic pulse is to introduce a reduction of the growth 
rates for the decay and modulational instabilities. Such spectral broadening could 
be due to a finite pulse coherence length, or through the use of random phase filters, 
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and would stabilize the propagation of electromagnetic pulses. 
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The physics of electron-positron plasmas is important for understanding the pulsar en- 
vironments and laboratory plasmas irradiated by intense lasers |3|, |4j , where quantum 



field effects come into play 5]. Pair plasmas are also believed to be important in the earl 
universe, in supernova remnants and active galactic nuclei, and in gamma-ray bursts 
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^ ■ A pair plasma can be created in these environments by collisions between strongly accel- 
erated particles 0. There is also a possibility of pair creation via high-energy curvature 
O | radiation photons in the vicinity of strongly magnetized astrophysical objects, triggered by 
charged particles streaming along the curved magnetic field 9|. This results in large quan- 



tities of positrons produced close to pulsar polar caps lfj,|lj|. In laboratory environments, 
laser experiments with focal spot intensities exceeding 10 20 W/cm 2 have demonstrated the 



production o 



interactions 



MeV electrons and evidence of positrons via laser multiphoton-gamma photon 
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I3I ] as well as via electron collisions 14, I5I ] 



When the intensity of the radiation is increased, the pair plasma particles attain rela- 
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ocities, and the plasma dynamics turns out to be essentially nonlinear, see 
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New decay instabilities, which are due to relativistic effects, can 
201. The nonlinear propagation of radiation in electron-positron plasmas has 



thus appear 

therefore been sugg ested to be the reason behind the high effective temperatures of pulsar 
radio emissions |21| as well as a possible source of the large frequency shifts in such emissions 

In the present paper, we will study the effects of partial coherence on the nonlinear 
radiation propagation in pair plasmas. It is shown that a spatial spectral broadening of the 
electromagnetic pulse can lead to a reduced growth rate for the decay and modulational 
instabilities. This type of pulse broadening could be due to a finite pulse coherence length 
(as in ultra-short laser pulse application), or through the use of random phase filters. The 
effect of the broadening is to stabilize the propagation of electromagnetic pulses in pair 
plasmas. 

The interaction between electromagnetic radiation and sound waves in pair plasmas with 
ultrarelativistic temperatures is given by the system of equations j^j 



(a t 2 -c 2 V 2 + 2^)A + ^( 



4 Tli 



3 n m 2 (T)c 4 



A 2 \A = 0, 



(1) 



and 



(d 2 - M 2 V 2 ) 



2^2 



ere 



:V 2 |A| 



(2) 



n 32T 2 

where A is the vector potential of the radiation, tu p = (7rnoe 2 c 2 /T) 1//2 is the effective plasma 
frequency, m(T) = 4T/c 2 is the effective mass, T is the electron and positron temperature, 
riQ is the unperturbed plasma density, ni (<^C n ) is the density perturbation, c is the speed 
of light in vacuum, e is the magnitude of the electron charge, u = [(c 2 /3)(l + aT 3 /Atiq)} 1 / 2 
is the sound speed in a relativistically hot pair plasma, a = An 2 /A5h 3 c 3 , and h is the Planck 
constant divided by 2n. A similar system has been previously derived for a cold plasma [24], 
and recently for a warm dust-laden plasma 

A kinetic description of the radiation propagation may be obtained by the introduction 
of the Wigner density matrix (see e.g. (26, 
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p ij (r,t,k,u) = -^ i J dTdtex P (ik-C-tLUT)(A*(r + C/2,t+T/2)A 3 (r-t/2,t-T/2)), (3) 



where the angular brackets denote the ensemble average (^J. The Wigner density matrix 
generalizes the the classical distribution function, and satisfies 



du dk p 
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where p = Tr(py) is the trace of the state density matrix. We note that for a monochromatic 
pump wave A = A exp(ik ■ r — iuj t) + c.c, we have / = 2| A \ 2 . As there is no polarization 
mixing and the ponderomotive force term only depends on J, it will here suffice to use p as 
our radiation field description. Applying the time derivative to the definition Q, using Eq. 
((TJ), and taking the trace of the resulting equation we obtain 



ujd t p + c 2 k ■ Vp + uj. 
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noting that the term proportional to the squared plasma frequency in Eq. (0) only contributes 
with a constant phase factor and therefore does not appear in Eq. (jSJ). The equations (J2J), 
(j3J), and (0) form a closed system describing the nonlinear interaction of sound waves and 
partially coherent radiation in pair plasmas with ultrarelativistic temperatures. 



Next, we analyze the modulational instability 0, 
(J2J), (JH), and (JHJ), by letting rii = hiexp(iK ■ r — iVLt) and p 
r — iQt), where \pi\ po an d K and fl are the perturbation wave vector and frequency, 



321 ] of the coupled system of equations 
p (fc,co>) + pexp(iK ■ 



respectively. We then linearize with respect to the perturbations to obtain the nonlinear 
dispersion relation 
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where po± = Po{k ± K/2, uj ± Q/2). For a monochromatic background distribution function 
p (k,u) = Io8(k — k )5(u — u> ), the expression (jHJ) reduces to 
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where 



D± = K 2 c 2 - tt 2 =f 2uj (tt -K-v g ), 
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and v g = k c 2 /u . We note that (J7J) agrees with (4) in Ref. |23j for decay processes when 
£>_ « 0. 

Next, we consider the case of partial coherence. We will here restrict our consideration 
to the one-dimensional case. Assuming that the background electromagnetic wave has a 
random phase (p(z) such that (exp[— i<p(z + C/2)] exy\iip(z — Q/2)]) = exp(— A|C|), the quasi- 



particle distribution takes the Breit-Wigner form 
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Po(fc '" ) = 7 (fc-fc ) 2 + A^ ( """ o) - (9) 
Here A is the half-width of the distribution centered around fco- Then the dispersion relation 

(0) can be written as 
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We next normalize to dimensionless variables according to uj$ — > ujo/uj p , Q — > Q/u> p , 
K — > Kc/up, A — > Ac/u p , u — > u/c, v g — > v g /c, and Jo — > Iqc 2 /24T 2 . In the case of the 
decay instability, we neglect the self-modulation of the electromagnetic wave, and we also 
have D_ w 0. Neglecting the fi 2 term in £)_, we have plotted the growth rate T of the 
decay instability in Fig. 1. In Fig. 2 we have plotted the growth rate for the modulational 
instability. In both Figs. 1 and 2 the damping effect of the spatial spectral broadening can 
be seen. We have chosen the normalized parameter values c^o = 1, lo = 0.1, v g = 0.1, and 
u = 0.2 in both Figs. 1 and 2. 

We also plot the case of an over-dense pair plasma, such that ujq = 0.2 in Fig. 3. The 
other parameters remain the same as above. 

In Fig. 4 we display the modulational instability growth rate for an under-dense pair 
plasma, such that ujq = 2. The remaining parameters remain the same as above. 

To summarize, we have consider the nonlinear propagation of spatial broadband elec- 
tromagnetic waves in a relativistically hot electron-positron plasma. For this purpose, we 
have introduced the Wigner transformation on a pair of equations comprising the nonlinear 
Schrodinger equation and a radiation pressure driven sound equation. As a result, we obtain 
a wave kinetic equation that is appropriate to derive a nonlinear dispersion relation. The 
latter exhibits that the growth rates of the decay and modulational instabilities are reduced 
when electromagnetic waves a broad spatial bandwidth. Hence, the electromagnetic waves 
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FIG. 1: The decay growth rate plotted as a function of K. The full curve has A = 0, while the 
dashed curve has A = 0.25. 



would be able to propagate over long distances without losing their energy. Our results 
should be useful in understanding the nonlinear propagation of partially coherent electro- 
magnetic waves in the early universe as well as in forthcoming laser-plasma interaction 
experiments. 
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FIG. 2: The modulational instability growth rate plotted as a function of K. The full curve has 
A = 0, while the dashed curve has A = 0.25. 
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FIG. 3: The modulational instability growth rate plotted as a function of K in the case of an 
over-dense pair plasma (loq = 0.2). The full curve has A = 0, while the dashed curve has A = 0.25. 
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FIG. 4: The modulational instability growth rate plotted as a function of K in the case of an 
under-dense pair plasma (ujq = 2). The full curve has A = 0, while the dashed curve has A = 0.25. 



